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Abstract. This is a tutorial for using a new g-series Maple pack- 
age. The package includes facilities for conversion between g-series 
and g-products and finding algebraic relations between g-series. 
Andrews found an algorithm for converting a g-series into a prod- 
uct. We provide an implementation. As an application we are 
able to effectively find finite q-product factorisations when they 
exist thus answering a question of Andrews. We provide other 
applications involving factorisations into theta-functions and eta- 
products. 
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1. Introduction 

In the study of g-series one is quite often interested in identifying 
generating functions as infinite products. The classic example is the 
Roger s-Ramanuj an identity 



DC 



q n -i-r 1 



(„. n\ II 



Here we have used the notation in (2.2). It can be shown that the 
left-side of this identity is the generating function for partitions whose 
parts differ by at least two. The identity is equivalent to saying the 
number of such partitions of n is equinumerous with partitions of n 
into parts congruent to ±1 (mod 5). 

The main goals of the package are to provide facility for handling 
the following problems. 

1. Conversion of a given g-series into an "infinite" product. 

2. Factorization of a given rational function into a finite g-product 
if one exists. 

3. Find algebraic relations (if they exist) among the g-series in a 
given list. 

A g-product has the form 

JV 

where the bj are integers. 

In 0, §10.7], George Andrews also considered Problems 1 and 2, 
and asked for an easily accessible implementation. We provide im- 
plementations as well as considering factorisations into theta-products 
and eta-products. The package provides some basic functions for com- 
puting g-series expansions of eta functions, theta functions, Gaussian 
polynomials and g-products. It also has a function for sifting out coeffi- 
cients of a g-series. It also has the basic infinite product identities: the 
triple product identity, the quintuple product identity and Winquist's 
identity. 



1.1. Installation instructions. The qseries package can be down- 
load via the WWW. First use your favorite browser to access the url: 
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http : / /www . math . uf 1 . edu/ ~f rank/qmaple . html then follow the di- 
rections on that page. There are two versions: one for UNIX and one 
for WINDOWS. 



2. Basic functions 

We describe the basic functions in the package which are used to 
build g-series. 

2.1. Finite g-products. 

2.1.1. Rising q-factorial. aqprod(a,q,n) returns the product 
(2.2) (a; q) n = (1 - a)(l - aq) ■ ■ ■ (1 - aq^ 1 ). 

We also use the notation 

oo 

(a-q) 00 = H(l-aq n - 1 ). 

n=l 

2.1.2. Gaussian polynomials. When < m < n, qbin(q,m,n) returns 
the Gaussian polynomial (or g-binomial coefficient) 



n 
m 



(q)n 



(<l)m(q), 



otherwise it returns 0. 
2.2. Infinite products. 

2.2.1. Dedekind eta products. Suppose > 0, and q = exp(2nir). 
The Dedekind eta function fZR p. 121] is defined by 



7]{t) = exp(7rzr/12) J^[(l — exp(27rmr)) 

n=l 

oo 

= q^\[{l-q n ). 

n=l 

etaq(q,k ,T) returns the g-series expansion (up to q T ) of the eta prod- 
uct 

oo 

l[(l-q kn )- 

n=l 

This corresponds to the eta function ^(fcr) except for a power of q. 
Eta products occur frequently in the study of g-series. For example, 
the generating function for p(n), the number of partitions of n, can be 
written as 

S p(n) '" = re .d-rr 
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See [I], pp. 3-4]. The generating function for the number of partitions 
of n that are p-cores [JIT] , a p (n), can be written 



E a * 

n=0 



n)q- 



n=l v 



pn\p 



q n ) 



Recently, Granville and Ono [|19|] were able to prove a long-standing 
conjecture in group representation theory using elementary and function- 
theoretic properties of the eta product above. 



2.2.2. Theta functions. Jacobi |22|, Vol I, pp. 497-538] defined four 



theta functions 6i(z,q) } i = 1,2,3,4. See also [^8[ Ch. XXI]. Each 
theta function can be written in terms of the others using a simple 
change of variables. For this reason, it is common to define 



9(z,q): = J2 *V- 

n=— oo 

theta(z,q,T) returns the truncated theta-series 

E *V- 

i=-T 

The case z = 1 of Jacobi's theta functions occurs quite frequently. We 
define 

oo 

fa(?) := E ^ +1)2/2 

n=— oo 
oo 

OM := E ^ 

n=— oo 
oo 

»4(?) := E (- 1 )"^ 2 

n=— oo 

theta2(q,t), theta3(q,t), theta4(q,t) (resp.) returns the g-series 
expansion to order O (q T ) of ^(g), 9z(q), 9i{q) respectively. Let a, and 
b be positive integers and suppose \q\ < 1. Infinite products of the form 



{q a ;q b )oo{q h - a ;q b U 



occur quite frequently in the theory of partitions and g-series. For ex- 
ample the right side of the Rogers- Ramanuj an identity is the reciprocal 
of the product with (a, b) = (1,5). The function jacprod(a,b,q,T) 
returns the g-series expansion to order 0(q T ) of the Jacob i-type infinite 
product given above. 
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3. Product Conversion 

In |], p. 233], || §10.7] there is a very nice and useful algorithm for 
converting a g-series into an infinite product. Any given g-series may 
be written formally as an infinite product 

oo oo 

i + X>^ n = II( 1 --<7T a "- 

n=l n=l 

Here we assume the b n are integers. By taking the logarithmic deriva- 
tive of both sides we can obtain the recurrence 

n 

nb n = b n -j da d . 
i =1 d\j 

Letting a n = 1 we obtain well-known special case 

n 

np(n) = ^2p(n-j)a(j). 

3=1 

We can also easily construct a recurrence for the a n from the recurrence 
above. 

The function prodmake is an implementation of Andrews' algorithm. 
Other related functions are etamake and jacprodmake. 

3.1. prodmake. prodmake (f ,q,T) converts the g-series / into an infi- 
nite product that agrees with / to 0(g T ). Let's take a look at the left 
side of the Roger's Ramanujan identity. 

> with(qseries) : 

> x:=l: 

> for n from 1 to 8 do 

x := x + q A (n*n)/aqprod(q,q,n) : 

od: 

> x := series(x,q,50) ; 

x ■= 1 + q + q 2 + q 3 + 2 g 4 + 2 q 5 + 3 q 6 + 3 q 7 + 4 q 8 + 5 q 9 + 6 g 10 
+7 g 11 + 9 g 12 + 10 g 13 + 12 g 14 + 14 g 15 + 17 g 16 + 19 g 17 + 23 g 18 
+26 g 19 + 31 g 20 + 35 g 21 + 41 g 22 + 46 g 23 + 54 g 24 + 61 g 25 + 70 
g 26 + 79 g 27 + 91 g 28 + 102 g 29 + 1 17 g 30 + 131 g 31 + 149 g 32 + 167 
g 33 + 189 g 34 + 211 g 35 + 239 g 36 + 266 g 37 + 299 g 38 + 333 g 39 + 
374 g 40 + 415 g 41 + 465 g 42 + 515 g 43 + 575 g 44 + 637 g 45 + 709 g 46 
+783 g 47 + 871 g 48 + 961 g 49 + 0(g 5 °) 

> prodmake (x,q, 40 ) ; 

1/ ((1 - g)(l - g 4 )(l - g 6 )(l - g 9 )(l - g n )(l - g 14 )(l - g 16 )(l - g 19 ) 
(1 - <? 21 )(1 " ? 24 )(1 - S 26 )(l - g 29 )(l - ^ 31 )(1 - ? 34 )(1 - ? 36 )(1 - ? 39 )) 
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We have rediscovered the right side of the Rogers- Ramanuj an identity! 
Exercise 1. Find (and prove) a product form for the g-series 



3.2. qfactor. The function qfactor is a version of prodmake. 
qf actor (f ,T) attempts to write a rational function / in q as a q- 
product, ie. as a product of terms of the form (1 — q l ). The second 
argument T is optional. It specifies an an upper bound for the expo- 
nents of q that occur in the product. If T is not specified it is given a 
default value of 4d + 3 where d is the maximum of the degree in q of the 
numerator and denominator. The algorithm is quite simple. First the 
function is factored as usual, and then it uses prodmake to do further 
factorisation into g-products. Thus even if only part of the function 
can be written as a g-product qfactor is able to find it. 

As an example we consider some rational functions T(r, h) intro- 
duced by Andrews [[|(p.l4) to explain Rogers's |32] first proof of the 



Rogers- 


Ramanuj an 


follows' 




(3.3) 


T(r,0)-- 


(3.4) 


T(r,l)-- 


(3.5) 


T(r, N) -- 



The T(r, n) are defined recursively as 



1, 
0. 



r + 2j 
3 



T(r + 2j,N-2j). 



> T:=proc(r,j) 

> option remember; 

> local x,k; 

> x:=0; 

> if j=0 or j=l then 

> RETURN((j-l) A 2) : 

> else 

> for k from 1 to floor(j/2) do 

> x:=x-qbin(q,k,r+2*k)*T(r+2*k, j-2*k) ; 

> od: 

> RETURN (expand (x) ) ; 

> fi: 

> end: 
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> t8:=T(8,8); 

t8 := 3 g 9 + 21 g 16 + 36 g 25 + 9 g 36 + g 6 + q 7 + 2 g 8 + 5 g 10 + 6 g 11 
+ 9 g 12 + 1 1 g 13 + 15 g 14 + 17 g 15 + 33 g 28 + 34 g 27 + 37 g 26 
+ 38 g 24 + 36 g 23 + 37 g 22 + 34 g 21 + 33 g 20 + 29 g 19 + 28 g 18 
+ 23 g 17 + 5 g 38 + 6 g 37 + 11 g 35 + 15 g 34 + 17 g 33 + 21 g 32 
+ 23g 31 + 28g 30 + 29 g 29 + 3g 39 + g 42 + g 41 + 2g 40 

> f actor (t8); 

g 6 (g 4 + g 3 + g 2 + g + l) (g 4 - g 3 + g 2 - g + l) 

(g 10 + g 9 + g 8 + g 7 + g 6 + g 5 + g 4 + g 3 + g 2 + g + l) (g 4 + 1) 
(g 6 + g 3 + 1) (g 8 + l) 

> qf actor (t8, 20) ; 

(l-g 9 )(l-g 10 )(l-g n )(l-g 16 )g 6 
(1 - g) (1 - g 2 ) (1 - g 3 ) (1 - g 4 ) 



Observe how we used factor to factor t8 into cyclotomic polynomials. 
However, qfactor was able to factor t8 as a g-product. We see that 

rp(o (g 9 ;g) 3 (i-g 16 )9 6 

1 (S, SJ = r . 



Exercise 2. Use qfactor to factorize T(r, n) for different 
values of r and n. Then write T(r, n) (defined above) as a 
g-product for general r and n. 

For our next example we examine the sum 

a + b 
b + k ' 

> dixson:=proc(a,b,c,q) 

> local x,k,y; 

> x:=0: y :=min(a,b,c) : 

> for k from -y to y do 

> x:=x+(-l) A k*q A (k*(3*k+l)/2)* 

> qbin(q,c+k,b+c)*qbin(q,a+k,c+a)*qbin(q,b+k,a+b) ; 

> od: 

> RETURN (x): 

> end: 

> dx := expand (dixson (5 , 5 , 5 , q) ) : 



(-l) k q k ( 3k+i y 2 



b + c 




c + a 


c + k 




a + k 
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> qf actor (dx, 20) ; 

(1 - g 6 ) (1 - g 7 ) (1 - g 8 ) (1 - g 9 ) (1 - g 10 ) (1 - g 11 ) (1 - g 12 ) 

(1 - g 13 ) (1 - g 14 ) (1 - g 15 )/((l - qf (1 - g 2 ) 2 (1 - g 3 ) 2 (1 - g 4 ) 2 

We find that 

' i<f'-- <l) : 



(3.6) 



V" (_l)*g*(3fc+l)/2 



k=— oo 



10 

5 + k 



110 



(95 9)1 ' 



Exercise 3. Write the sum 



fc(3fc+l)/2 



9 

k=— oo 

as a g-product for general integral a 



2a 
a + k 



3.3. etamake. Recall from ( |2.2.1| ) that etaq is the function to use for 
computing g-expansions of eta products. If one wants to apply the the- 
ory of modular forms to g-series it is quite useful to determine whether 
a given g-series is a product of eta functions. The function in the pack- 
age for doing this conversion is etamake. etamake (f ,q,T) will write 
the given g-series / as a product of eta functions which agrees with / 
up to q T . As an example, let's see how we can write the theta functions 
as eta products. 

> theta2(q,100)/q A (l/4) ; 

2g 132 + 2g 110 + 2g 90 + 2g 72 + 2g 56 + 2g 42 + 2g 30 + 2g 20 + 2g 12 + 2g 6 
+ 2 g 2 + 2 + g 156 

> etamake (" ,q, 100) ; 

o ^(4r) 2 
g 1 /4 r? (2r) 

> theta3(q,100) ; 

2g 121 + 2g 100 + 2g 81 + 2g 64 + 2g 49 + 2g 36 + 2g 25 + 2g 16 + 2g 9 + 2g 4 
+ 2g+ 1 

> etamake (" ,q, 100) ; 

r](Ar) 2 i]( t ) 2 
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> theta4(q,100) ; 

-2g 121 + 2g 100 -2g 81 + 2g 64 -2g 49 + 2g 36 -2g 25 + 2g 16 -2g 9 + 2g 4 
-2q+l 

> etamakeC ,q, 100) ; 

V(r) 2 
rj(2r) 

We are led to the well-known identities: 

^(4r) 2 



e 2 (q) = 2 



r,(2r) ' 
r?(2r) 5 
?7(4r) 2 ?7(r) 2 ' 



r)(2ry 

The idea of the algorithm is quite simple. Given a g-series / (say with 
leading coefficient 1) one just keeps recursively multiplying by powers 
of the right eta function until the desired terms agree. For example, 
suppose we are given a g-series 

/ = 1 + cq k + • • • . 

Then the next step is to multiply by etaq(q,k,T) A (-c). 



Exercise 4. Define the g-series 



oo oo 

q" 



a{q) ■= ' 



n=— oo m=— oo 

oo oo 



b(q) := E E 



q 



n=—oo m=—oo 
oo oo 



c ( 5 ) ;= - (ri+i/3)2+( " +i/3)(m+i/3)+{m+i/3i:: 



q y 

n=— oo m=—oo 



where u = exp(27ri/3). Two of the three functions above can 
be written as eta products. Can you find them? 
Hint: It would be wise to define 

> omega := Root0f(z A 2 + z + 1 = 0) ; 
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3.4. jacprodmake. In (|2.2.2| ) we observed that the right side of the 
Roger s-Ramanuj an identity could be a written in terms of a Jacobi 
product. The function jacprodmake converts a g-series into a Jacobi- 
type product if one exists. Given a g-series /, jacprodmake (f ,q,T) 
attempts to convert / into a product of theta functions that agrees 
with / to order 0(g T ). Each theta-function has the form JAC(a, b, oo), 
where a, b are integers and < a < b. If < a, then JAC(a, b, oo) 
corresponds to the theta-product 

(q a ;q b U(q b ~ a ;q b )oo(q b ;q b )oo. 

We call this a theta product because it is #(— g^ 2 ^/ 2 , g b//2 ). The 
jacprodmake function is really a variant of prodmake. It involves using 
prodmake to compute the sequence of exponents and then searching for 
periodicity. 

If a = 0, then JAC(0, b, oo) corresponds to the eta-product 

(?W)oc 

Let's re-examine the Roger s-Ramanuj an identity. 

> with(qseries) : 

> x:=l: 

> for n from 1 to 8 do 

> x : =x+q A (n*n) / aqprod (q , q , n) : 

> od: 

> x:=series(x,q,50) : 

> y :=jacprodmake(x,q,40) ; 

X4C(0,5,oo) 
V ''~ X4C(l,5,oo) 

> z:=jac2prod(y) ; 

1 

Z '. = 

0,9 5 )oo(g 4 ,g 5 )oo 

Note that we were able to observe that the left side of the Rogers- 
Ramanujan identity (at least up through g 40 ) can be written as a 
quotient of theta functions. We used the function jac2prod, to sim- 
plify the result and get it into a more recognizable form. The func- 
tion jac2prod(jacexpr) converts a product of theta functions into 
g-product form; ie. as a product of functions of the form (g a ;g 6 )oo- 
Here jacexpr is a product (or quotient) of terms JAC(i,j, oo), where 
i, j are integers and < i < j. 

A related function is jac2series. This converts a Jacobi-type prod- 
uct into a form better for computing its g-series. It simply replaces each 
Jacobi-type product with its corresponding theta-series. 



> with(qseries) 
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> x:=0: 

> for n from to 10 do 

x := x + q A (n*(n+l)/2)*aqprod(-q,q,n)/aqprod(q,q,2*n+l) 
od: 

> x:=series(x,q,50) : 

> jacprodmake(x,q,50) ; 

JAC(0, 14, oo) 6 /(JAC(1, 14, oo) 2 JAC(3, 14, oo) JAC(4, 14, oo) 



JAC(5,14,oo,JAC(6, 14, oo) ^Og 14, g } 
> jac2series(" ,500) ; 

( g 364 _ g 210 + g 98 _ g 28 + 1 _ g 14 + g 70 _ g 168 + g 308 _ g 490 )6 / (( 

_ g 621 + g 496 _ g 385 + g 288 _ g 205 + g 136 _ g 81 + g 40 _ g 13 + j _ g 
+ gl6 _ g 45 + g 88 _ g 145 + g 216 _ g 301 + g 400 _ g 513 )2( _ g 603 + g 480 
_ g 371 + g 276 _ g 195 + g 128 _ g 75 + g 36 _ g ll + x _ g 3 + g 20 _ g 51 
+ g 96 _ g 155 + g 228 _ g 315 + g 416 _ g 531 )( _ g 594 + g 472 _ g 364 + g 270 
_ g 190 + g 124 _ g 72 + g 34 _ g W + 1 _ g 4 + g 22 _ g 54 + g 100 _ g 160 
+ g 234 _ g 322 + g 424 _ g 540 )( _ g 585 + g 464 _ g 357 + g 264 _ g 185 + g 120 
_ g 69 + g 32 _ g 9 + 1 _ g5 + g 24 _ g 57 + g 104 _ g 165 + g 240 _ g 329 

_|_ g 432 _ g 549 )(_g 576 _|_ g 456 _ g 350 _|_ g 258 _ g 180 _|_ g 116 _ g 66 _|_ g 30 

_ g 8 + 1 _ g6 + g 2 6 _ g 60 + g 108 _ g 170 + g 246 _ g 336 + g 440 _ g 558 )(( 

-2 g 567 + 2 g 448 - 2 g 343 + 2 g 252 - 2 g 175 + 2 g 112 - 2 g 63 + 2 g 28 
-2g 7 +l)/( 

g 3 6 4 _ g 2 10 + g 98 _ g 28 + 1 _ gl4 + g 70 _ g 168 + g 308 _ g490 )) l/2 ) 

It seems that the g-series 



(~g; q) n q n{n+1)/2 



n>0 



can be written as Jacob i-type product. Assuming that this is the case 
we used jac2series to write this g-series in terms of theta-series at 
least up to g 500 . This should provide an efficient method for computing 
the g-series expansion and also for computing the function at particular 
values of q. 



Exercise 5. Use j acprodmake and jac2series to compute 
the g-series expansion of 

-g; g)„g™( 3n+1 )/ 2 
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up to g 1000 , assuming it is Jacobi-type product. Can you iden- 
tify the infinite product? 



4. The Search for Relations 

The functions for finding relations between g-series are findhom, 
f indhomcombo, f indnonhom, f indnonhomcombo, and f indpoly. 

4.1. findhom. If the g-series one is concerned with are modular forms 
of a particular weight, then theoretically these functions will satisfy 
homogeneous polynomial relations. See |16|, p. 263], [^] for more details 
and examples. The function f indhom(L,q,n,topshif t) returns a set 
of potential homogeneous relations of order n among the g-series in the 
list L. The value of topshift is usually taken to be zero. However 
if it appears that spurious relations are being generated then a higher 
value of topshift should be taken. 

The idea is to convert this into a linear algebra problem. This pro- 
gram generates a list of monomials of degree n of the functions in the 
given list of g-series L. The g-expansion (up to a certain point) of each 
monomial is found and converted into a row vector of a matrix. The 
set of relations is then found by computing the kernel of the transpose 
of this matrix. As an example, we now consider relations between the 
theta functions 9 3 (q), ^(g), # 3 (g 2 ), and # 4 (g 2 ). 

> with(qseries) : 

> findhom([theta3(q,100) ,theta4(q, 100) ,theta3(q A 2 , 100) , 
theta4(q A 2,100)] ,q,l,0) ; 

# of terms , 25 

RELATIONS of order—, 1 

{{}} 

> findhom([theta3(q,100) ,theta4(q, 100) ,theta3(q A 2 , 100) , 
theta4(q A 2,100)] ,q,2,0) ; 

# of terms , 31 

RELATIONS of order—, 2 

{X, 2 + X 2 2 - 2X 3 2 , -X ± X 2 + X 4 2 } 

From the session above we see that there is no linear relation between 
the functions # 3 (g), 6* 4 (g), 6* 3 (g 2 ) and # 4 (g 2 ). However, it appears that 
there are two quadratic relations: 

and 



2 

1/2 



Hf) = (ol{ q )el{ q )y\ 
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This is Gauss' parametrization of the arithmetic-geometric mean iter- 
ation. See Ch. 2] for details. 



Exercise 6. Define a(q), b(q) and c(q) as in Exercise 2. 
Find homogeneous relations between the functions a(q), b(q), 
c(q), a(q 3 ), b(q 3 ), c(q 3 ). In particular, try to get a(q 3 ) and 
b(q 3 ) in terms of a(q), b(q). See |I(J for more details. These 
results lead to a cubic analog of the AGM due to Jon and Peter 
Borwein §, @. 



4.2. f indhomcombo. The function f indhomcombo is a variant of f indhom. 
Suppose / is a g-series and L is a list of g-series. f indhomcombo (f , L , q, 
n,topshif t ,etaoption) tries to express / as a homogeneous polyno- 
mial in the members of L. If etaoption=yes then each monomial in 
the combination is "converted" into an eta-product using etamake. 

We illustrate this function with certain Eisenstein series. For p an 
odd prime define 

xfa) — I — ) (the Legendre symbol). 



V 

Suppose k is an integer, k > 2, and (p— l)/2 = k (mod 2). Define the 
Eisenstein series 



m=l n=l 



Then U Pt k is a modular form of weight k and character \ for the con- 
gruence subgroup r (p). See |f?6|j , [0 for more details. The classical 



result is the following identity found by Ramanujan |30], Eq. (1.52), p. 
354]: 

?7(5r) 5 



5,2 



7?(r) 



Kolberg [2^j has found many relations between such Eisenstein series 
and certain eta products. The eta function r/(r) is a modular form of 
weight | p5| , p. 121]. Hence the modular forms 



B = r/(5r) 5 B = y{r) 5 
1 ' ' 2 ' ?](5r) 

are modular forms of weig ht = 2. In fact, it can be shown that 
they are modular forms on To (5) with character fgj. We might there- 
fore expect that U^^ can be written as a homogeneous cubic polynomial 
in Bi, B2. We write a short maple program to compute the Eisenstein 
series U p> k- 
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> with(numtheory) : 

> UE:=proc(q,k,p, trunk) 

> local x,m,n: 

> x : =0 : 

> for m from 1 to trunk do 

> for n from 1 to trunk/m do 

> x:=x + L(m,p)*n A (k-l)*q A (m*n) : 

> od: 

> od: 

> end: 

The function UE(q,k,p, trunk) returns the g-expansion of U Pt k up 
through gtrunk ^y e no ^ e ^at L(m,p) returns the Legendre symbol 
m\ 

— . We are now ready to study U5 e- 
V ) 

> with(qseries) : 

> f := UE(q, 6,5,50) : 

> Bl := etaq(q,l,50) A 5/etaq(q,5,50) : 

> B2 := q*etaq(q,5,50) A 5/etaq(q,l,50) : 

> f indhomcombo(f , [B1,B2] ,q, 3,0, yes) ; 

# of terms , 25 
possible linear combinations of degree , 3 

r](5 r) 3 r/(r) 9 + 40 r/(5 t) 9 t](t) 3 + 335 

7]{tY 

{X l 2 X 2 + 40XxX 2 2 + 335 X 2 3 } 
It would appear that 

^5,6 = r/(5 r f^r f + 40 r/(5 r) 9 r ? (r) 3 + 335 

The proof is a straightforward exercise using the theory of modular 
forms. 

Exercise 7. Define the following eta products: 

rj{T) 7]{7t) 
What is the weight of these modular forms? 

Write Ur t 3 in terms of C±, C2, C3. 

The identity that you should find was originally due to Ra- 
manujan. Also see Fine |13|, p. 159] and Eq. (5.4)]. 
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4.3. f indnonhom. In section f4.1| we introduced the function findhom 
to find homogeneous relations between g-series. The nonhomogeneous 
analog is f indnonhom. 

The syntax of f indnonhom is the same as findhom. Typically (but 
not necessarily) findhom is used to find relations between modular 
forms of a certain weight. To find relations between modular functions 
we would use f indnonhom. We consider an example involving theta 
functions. 

> with(qseries) : 

> F := q -> theta3(q,500)/theta3(q A 5,100) : 

> U := 2*q*theta(q A 10,q A 25,5)/theta3(q A 25,20) ; 

2 q (g 575 + g 360 + q 195 + q 80 + q 15 + 1 + g 35 + q 120 + g 255 + g 440 + q 675 ) 



U : 



2 g 625 + 2 g 400 + 2 g 225 + 2 g 100 + 2 g 25 + 1 

> EQNS := f indnonhom ([F(q) ,F(q A 5) ,U] ,q,3,20) ; 

# of terms , 61 

RELATIONS of order—, 3 

EQNS := {-1 - X!X 2 X 3 + X 2 2 + X 3 2 + X 3 } 

> ANS:=EQNS[1] ; 

ANS := -1 - X!X 2 X 3 + X 2 2 + X 3 2 + X 3 

> CHECK := subs(X[l]=F(q) ,X[2]=F(q A 5) ,X[3]=U,ANS) 

> series (CHECK, q, 500) ; 

O (g 500 ) 



We define 

F(,) r " M 



and 



U(q) :=2 



00 

25n 2 +10n+l 



. n=— 00 



3 (q 25 ) 

We note that U (g) and F(q) are modular functions since they are ratios 
of theta series. From the session above we see that it appears that 

1 + F(q)F(q 5 )U(q) = F(q 5 f + U(q) 2 + U(q). 
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Observe how we were able to verify this equation to high order. When 
f indnonhom returns a set of relations the variable X has been declared 
global. This is so we can manipulate the relations. It this way we were 
able to assign ANS to the relation found and then use subs and series 
to check it to order (g 500 ). 

4.4. f indnonhomcombo. The syntax of f indnonhomcombo is the same 
as f indhomcombo. We consider an example involving eta functions. 
First we define 

t;(49 r) 
rj(r) 

and 




V(r) J 



Using the theory of modular functions it can be shown that one must 
be able to write T 2 in terms of T and powers of £. We now use 
f indnonhomcombo to get T 2 in terms of £ and T. 

> xi:=q A 2*etaq(q, 49, 100)/etaq(q, 1,100) : 

> T:=q*(etaq(q,7, 100) /etaq(q, 1 , 100) ) A 4 : 

> f indnonhomcombo (T A 2, [T,xi] ,q,7,-15) ; 

# of terms , 42 
matrix is , 37, x, 42 
possible linear combinations of degree , 7 

{ 147 X 2 5 + 343 X 2 7 + 343 X 2 6 + X 2 + 49 X 2 4 + 49 X X X 2 Z + 7 X X X 2 
+21 X 2 3 + 7X 2 2 + 35X 1 X 2 2 } 



Then it seems that 

T 2 = (35 i 2 + 49 £ 3 + 7f) T+343£ 7 +343£ 6 +147£ 5 +49£ 4 +21£ 3 +7£ 2 +£. 

This is the modular equation used by Watson [|3~7|1 to prove Ramanujan's 
partition congruences for powers of 7. Also see H and p4|, and see 



1% for an elementary treatment. 
Exercise 8. Define 



»j(25t) 



and 

t](5t 



T :-- 



r}(r) 

Use f indnonhomcombo to express T as a polynomial in £ of 
degree 5. The modular equation you find was used by Watson 
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to prove Ramanujan's partition congruences for powers of 5. 
See [21]] for an elementary treatment. 



Exercise 9. Define a(q) and c(q) as in Exercise 2. Define 

/ \ c (q) 3 

a(q) 



and the classical Eisenstein series (usually called E 3 ; see |J3], p. 
93]) 



N(q) := 1 - 504 ^ T 



n°q 



Q 



n=l 

Use f indnonhomcombo to express N(q) in terms of a(q) and 
ac(g). 

HINT: N(q) is a modular form of weight 6 and a(q) and c(q) 
are modular forms of weight 1. See for this result and many 
others. 



4.5. f indpoly. The function f indpoly is used to find a polynomial 
relation between two given g-series with degrees specified. 

f indpoly(x,y,q,degl,deg2, check) returns a possible polynomial 
in X, Y (with corresponding degrees degl, deg2) which is satisfied by 
the g-series x and y. 

If check is assigned then the relation is checked to O ^g c heck J 

We illustrate this function with an example involving theta functions 
and the function a(q) and c(q) encountered in Exercises 2 and 7. It 
can be shown that 



a(q) = 9 3 (q)6 3 (q 3 ) + 9 2 (q)9 2 (q :i ). 



See |T(J for details. This equation provides a better way of computing 
the g-series expansion of a(q) than the definition. In Exercise 2 you 
would have found that 

c{q) = 3 . 

r/(r) 



See |TD] for a proof. Define 



c 3 



and 
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We use f indpoly to find a polynomial relation between x and y. 

> with(qseries) : 

> xl := radsimp(theta2(q,100) A 2/theta2(q A 3,40) A 2) : 

> x2 := theta3(q,100) A 2/theta3(q A 3,40) A 2: 

> x := xl+x2: 

> c := q*etaq(q,3,100) A 9/etaq(q,l,100) A 3: 

> a := radsimp(theta3(q,100)*theta3(q A 3,40)+theta2(q,100) 
*theta2(q A 3,40)) : 

> c := 3*q A (l/3)*etaq(q,3,100) A 3/etaq(q,l,l00) : 

> y := radsimp(c A 3/a A 3) : 

> Pl:=findpoly(x,y,q,3,l,60) ; 
WARNING: X,Y are global. 



dims , 8, 18 
The polynomial is 



(X + 6) 3 F-27 (X + 2) 2 



Checking to order, 60 



o (g 59 ) 

PI := (X + 6) 3 F-27 (X + 2) 2 



It seems that x and y satisfy the equation 

p(x, y) = {x + 6) 3 y - 27{x + 2) 2 = 0. 



Therefore it would seem that 



c 3 _ 2? (x + 2) 2 
a 3 (x + 6) 3 ' 
See ||, pp. 4237-4240] for more details. 



Exercise 10. 



Define 




Use polyf ind to find i/ = ^asa rational function in m. 



5. Sifting coefficients 



Suppose we are given a g-series 



oo 



n=0 
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Occasionally it will turn out the generating function 

oo 
n=0 

will have a very nice form. A famous example for p(n) is due to Ra- 
manujan: 



1 - g 5n ) 6 

n\5 ' 



n=0 n=l ^ y ' 

See 10, Cor. 10.6]. In fact, G.H. Hardy and Major MacMahon g!§ p. 
xxxv] both agreed that this is Ramanujan's most beautiful identity. 
Suppose s is the g-series 

^a^ + 0(g T ) 

then sift (s,q,n,k,T) returns the g-series 

Y^ a ^kq l + 0{q T l n ). 

We illustrate this function with another example from the theory of 
partitions. Let pd(n) denote the number of partitions of n into distinct 
parts. Then it is well known that 



oo oo oo , 

EM»)?" = II(i+? B ) = Il7r 



q 2n ) 



q n ) 

n=0 n=l ra=l v ^ ' 

We now examine the generating function of pd(5n + 1) in maple. 

> PD:=series(etaq(q,2,200)/etaq(q,l,200) ,q,200) : 

> sift(PD,q,5, 1,199) ; 

1 + 4 q + 5010688 g 26 + 53250 g 15 + 668 g 7 + 12 g 2 + 165 g 5 
+ 12076 g 12 + 1087744 g 22 + 109420549 g 35 + 76 g 4 + 32 g 3 
+ 340 g 6 + 1260 g 8 + 2304 g 9 + 4097 g 10 + 7108 g 11 + 20132 g 13 
+ 32992 g 14 + 84756 g 16 + 133184 g 17 + 206848 g 18 + 317788 g 19 
+ 728260 g 21 + 20792120 g 30 + 2368800 g 24 + 483330 g 20 
+ 1611388 g 23 + 3457027 g 25 + 7215644 g 27 + 10327156 g 28 
+ 14694244 g 29 + 29264960 g 31 + 40982540 g 32 + 57114844 g 33 
+ 79229676 g 34 + 150473568 g 36 + 206084096 g 37 
+ 281138048 g 38 + 382075868 g 39 

> PD1: = " : 

> etamake(PDl,q,38) ; 

r/(5r) 3 r/(2r) 2 
g 5 / 24 r]( 10 r)r](r) 4 
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So it would seem that 



OO OO 

J2pd(5n + l)q n = H [ —^ 

n=0 n=l ^ 



Y 5n\3/ 



1 — q 



2n\2 



q 10n )(l-q 



n\4 



This result was found originally by R0dseth [31 



Exercise 11. Rodseth also found the generating functions 
for pd(5n + r) for r = 0, 1, 2, 3 and 4. For each r use sift and 
jacprodmake to identify these generating functions as infinite 
products. 



6. Product Identities 



At present, the package contains the Triple Product identity, the 
Quintuple Product identity and Winquist's identity. These are the 
most commonly used of the Macdonald identities |28[ 



34 , 35 . The 



Macdonald identities are the analogs of the Weyl denominator for affine 
roots systems. Hopefully, a later version of this package will include 
these more general identities. 



6.1. The Triple Product Identity. The triple product identity is 

oo oo 

(6.7) (-i)vv (n - i)/2 = n^ 1 - - <? n /^)(i - 

n=— oo ?i=l 

where z ^ and \q\ < 1. The Triple Product Identity is originally 
due to Jacobi [^2], Vol I]. Andrews [0 and Lewis [^7| have found nice 
combinatorial proofs. The triple product occurs frequently in the the- 
ory of partitions. For instance, most proofs of the Roger s-Ramanuj an 
identity crucially depend on the triple product identity. 

tripleprod(z,q,T) returns the g-series expansion to order 0(q T ) 
of Jacobi's triple product ( |6.7|) . This expansion is found by simply 
truncating the right side of (|6.7p. 

> tripleprod(z,q, 10) ; 

+ ^--^ + ^-?- + l-z + z 2 q-zy + zy-z 5 q 10 + z (i q 15 
z° z° z^ z a z z z 

> tripleprod(q,q A 3, 10) ; 

q 57 - q 40 + q 26 - q 15 + q 7 - q 2 + 1 - q + q 5 - q 12 + q 22 - q 35 + q 51 
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The last calculation is an illustration of Euler's Pentagonal Number 
Theorem g], Cor. 1.7 p.ll]: 

(6.8) 

oo oo oo 

- f) = - g3n-l )(1 _ g 3n-2 )(1 _ ? 3 n) = ^ (^n^ (3n-l)/2_ 

n=l n=l ra=— oo 

6.2. The Quintuple Product Identity. The following identity is 
the Quintuple Product Identity: 

(6.9) (-z, q)oo(-~, q)oo(z 2 q, g 2 )oo(4, <? 2 )oo(<?, q)oo 



z z^ 



En \—3m / \3m-l\ 



3m-l\ _ m ( 3m+1 ) 



m=— oo 



Here | lg| < 1 and 2^0. This identity is the BC\ case of the Macdonald 
identities |28| . The quintuple product identity was is usually attributed 
to Watson []36|. However it can be found in Ramanujan's lost notebook 
[ p0| , p. 207]. Also see || for more history and some proofs. 

The function quinprod(z,q,T) returns the quintuple product iden- 
tity in different forms: 

(i) If T is a positive integer it returns the g-expansion of the right 
side of ( |6.9| ) to order 0(q T ). 

(ii) If T = prodid then quinprod(z,q, prodid) returns the quintuple 
product identity in product form. 

(hi) If T = seriesid then quinprod(z,q, seriesid) returns the quin- 
tuple product identity in series form. 

> quinprod (z , q , prodid) ; 



{-Z, q)oo ( , g)oo (z q, q )oo ( — , q ) 0O (q, q) c 

z z z 



2 

%,q 3 )oo (qz 3 ,q 3 )oo (q 3 ,q 3 )oo + z (\, q 3 )^ {q 2 z 3 ,q 3 ) OQ (q 3 ,q 3 ) 
z 6 z 6 



> quinprod(z, q, seriesid) ; 



-z, q)oo (--, q)oc {z 2 q,q 2 ) 00 q 2 )oo {q, q)c 

z z z 

oo 

v(3m-l)\ (l/2m(3m+l)) 



((-^)(- 3m )-(-^)( 3m - 1 ))g( 



> quinprod(z,q,3) ; 
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Let's examine a more interesting application. Euler's infinite product 
may be dissected according to the residue of the exponent of q mod 5: 

oo 

- q n ) = E (q) + qE^q 5 ) + q 2 E 2 (q 5 ) + q 3 E 3 (q 5 ) + q 4 E 2 (q 5 ). 

n=l 

By ( |6.8[ ) we see that E% = E4 = since n(3n — l)/2 = 0, 1 or 2 mod 
5. Let's see if we can identify E . 

> with(qseries) : 

> EULER:=etaq(q, 1,500) : 

> E0:=sift(EULER,q,5,0,499) ; 

E0:=l + q-q 3 -q 7 -q s -q l4 + q 20 + q 29 + q 31 + q 42 - q 52 - q m 
_ g 69_ g 85 + g 99 

> jacprodmake(E0,q,50) ; 

JAC(2,5,oo) JAC(0,5,oo) 



JAC(l,5,oo) 

> jac2prod("); 

(g 5 ,g 5 )oo (q 2 ,q 5 )oo (g 3 ,9 5 )oo 

> quinprod(q,q A 5,20) : 

> series (" ,q, 100) ; 

i+q-q 3 -q 7 - q 8 -q u +q 20 +q 29 +q 31 +q 42 -q 52 -q m -q m -q 85 +o{q 99 ) 

From our maple session it appears that 
(6 10) E = (^ 5 )°°(g 2 'g 5 )°°(g 3 ^ 5 )°° 

(q;q 5 )oo(q 4 ;q 5 )oc 

and that this product can be gotten by replacing q by q 5 and z by q in 
the product side of the quintuple product identity (|6.9|) . 



Exercise 12. (i) Use the quintuple product identity fl6.9|) 
and Euler's pentagonal number theorem to prove Q6.10D above. 

(ii) Use maple to identify and prove product expressions for 
Ei and E 2 . 

(iii) This time see if you can repeat (i), (ii) but split the 
exponent mod 7. 

(iv) Can you generalize these results to arbitrary modulus? 
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6.3. Winquist's Identity. Back in 1969, Lasse Winquist |39| discov- 
ered a remarkable identity 

(6.11) 

(a; q)oo(q/a; q)oo(b; q)oo(q/b; q)oo(ab; q)oo(q/ (a&); g)oo(a/&; g)oo(ty («<?); <?)oo(<?; q) 

oo oo 

= £ J] (_l)^(( a -3«_ a 3n+3 )(r 3m_ 6 3m+l ) 
n=0 m=— oo 

_|_ / a ~3m+l _ Q 3m+2w^3n+2 _ ^-3n-l\\^3n(n+l)/2+m(3m+l)/2 

By dividing both sides by (1 — a)(l — 6) and letting a, 6 — > 1 he was 

oo 

able to express the product J^[(l — q n ) 10 as a double series and prove 

n=l 

Ramanujan's partition congruence 

p(lln + 6) = (mod 11). 

This was probably the first truly elementary proof of Ramanujan's 
congruence modulo 11. The interested reader should see Dyson's article 



E| for some fascinating history on identities for powers of the Dedekind 



eta function and how they led to the Macdonald identities. A new proof 
of Winquist's identity has been found recently by S-Y Kang Mike 
Hirschhorn |20| has found a four-parameter generalization of Winquist's 
identity. 

The function winquist(a,b,q,T) returns the g-expansion of the 
right side of (|6.11| ) to order 0(q T ). 

We close with an example. For 1 < k < 33 define 

oo 

g(*) = n(l-g*)(l-g JB -*)(l-? 38 )- 

n=l 

Now define the following functions: 

A = Q(15), A3 = Q(12), A 7 = Q(Q), A 8 = Q(3), A 9 = Q(9); 

B = Q(16)-q 2 Q(5), 
B 1 = Q(U)-qQ(8), 
B 2 = Q(13)-q 3 Q(2), 
B 4 = Q(7) + qQ(A), 
B 7 = Q(10) + q 3 Q(l). 

These functions occur in Theorem 6.7 of ]15| as well as the function 
A B 2 - q 2 A 9 B 4 . 

> with(qseries) : 

> Q:=n->tripleprod(q A n,q A 33,10) : 

> A0:=Q(15): A3:=Q(12): A7:=Q(6): A8:=Q(3): A9:=Q(9): 
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> B2:=Q(13)-q A 3*Q(2) : B4:=Q(7)+q*Q(4) : 

> IDG:=series( ( A0*B2-q A 2*A9*B4) ,q,200) : 

> series(IDG,q, 10) ; 

1 - g 2 - 2 g 3 + g 5 + g 7 + g 9 + O (q 11 ) 

> jacprodmake(IDG,q,50) ; 

JAC( 2, 11, oo ) JAC( 3, 11, oo f JAC( 5, 11, oo ) 
JAC(0,ll,oo) 3 

> jac2prod("); 

(g 2 ,g n )oo (g 9 ,g n )oo (g n ,g u )oo (g 3 ,g n )oo 2 (g s ,q u )J (g 5 ,g n )oo 
(sV 1 )*, 

> series (winquist(q A 5,q A 3,q A 11, 10) ,q,20) ; 

l-g 2 -2g 3 + g 5 + g 7 + 0(g 9 ) 

> series(IDG-winquist(q A 5,q A 3,q A ll,10) ,q,60) ; 

0(g 49 ) 

From our maple session it seems that 
(6.12) 

A B 2 - q 2 A g B 4: = (g 2 ; g^Ug 9 ; q^)^ q U )M; 1 U Ul 8 ; ^)lo 

and that this product appears in Winquist's identity on replacing q by 
g 11 and letting a = g 5 and b = g 3 . 

Exercise 13. (i) Prove (|6.12|) by using the triple prod- 
uct identity (|6.7f) to write the right side of Winquist's identity 
( |6.11| ) as a sum of two products. 

(ii) In a similar manner find and prove a product form for 

A B - q 3 A 7 B^. 
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